STABILITY ESTIMATES FOR AN INVERSE SCATTERING 
PROBLEM AT HIGH FREQUENCIES 



HABIB AMMARI, HAJER BAHOURI, DAVID DOS SANTOS FERREIRA, 
AND ISABELLE GALLAGHER 



Abstract. We consider an inverse scattering problem and its near-field 
approximation at high frequencies. We first prove, for both problems, 
Lipschitz stability results for determining the low-frequency component 
of the potential. Then we show that, in the case of a radial poten- 
tial supported sufficiently near the boundary, infinite resolution can be 
achieved from measurements of the near-field operator in the monotone 
case. 
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Introduction 

The first aim of this paper is to establish Lipschitz stability results for 
the inverse scattering problem of determining the low-frequency component 
(lower than the operating frequency) of the compactly supported potential 
from scattering or near-field measurements. It is known that, in general, 
the problem is exponentially unstable [H El EU HQ]. However, taking 
advantage of a priori information may improve stability and give accurate 
reconstruction algorithms [8[ I12[, fT3] . The Lipschitz stability results proved 
in this paper together with the recent analysis of the local convergence of 
the nonlinear Landweber iteration in [33] show that the low-frequency com- 
ponent of the potential can be determined from the data in a linearly stable 
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way. Moreover, they precisely quantify the resolution limit, which is de- 
fined as the characteristic size of the smallest oscillations in the potential 
that can be stably recovered from the data. Since Rayleigh's work, it has 
been admitted that the resolution limit in inverse scattering is of order it 
over the operating frequency [19J. This is nothing else than a direct appli- 
cation of the uncertainty principle in inverse scattering |T8l Ell EEl [55] . It 
is well-known that if the support of the potential is a point support, then 
the reconstructed location of the point potential from the scattering data 
has finite size of order of the Rayleigh resolution limit [18] . Having this 
in mind, the results of this paper prove that the Fourier transform of the 
potential can be reconstructed in a linearly stable way for all frequencies 
(dual variable to the space one) smaller than the operating frequency, and 
therefore, justify the notion of resolution limit. More intriguingly, again 
in view of [33], they prove that the stability of the reconstruction of the 
potential increases at high operating frequencies. 

The second aim of the paper is to show that infinite resolution can be 
achieved from near-field measurements. Here, the near-field operator ap- 
proximates Sommerfeld's radiation condition and is equivalent to the mea- 
surements of the Cauchy data at a finite distance. Moreover, if the potential 
is supported near the boundary, then infinite resolution can be achieved in 
the monotone case. In fact, a Lipschitz stability result holds for both the 
low and high frequency components of the potential. It should be noted 
that the scattering amplitude can be recovered from the near-field operator. 
However, approximating the near-field operator from the scattering ampli- 
tude is a severely ill-posed problem [371 S3 HS1 ES] and therefore cannot 
be of any practical and realistic use. It was shown in [37] that in order to 
compute the near-field operator from the scattering amplitude one needs to 
differentiate the scattering amplitude an infinite number of times. 

The results of this paper extend to medium scattering the recent results in [91 
[TOl [63] , where a stability and resolution analysis was performed for linearized 
conductivity and wave imaging problems. They can be also used to justify 
the hopping (or continuation in the frequency) reconstruction algorithms 
proposed in [11 [Ml ESI [26] . 

In connection with our results in this paper, we also refer to the works by 
Isakov [3l] and Isakov and Kindermann [35|, Bao, Lin, and Triki [151 116j . 
Nagayasu, Uhlmann and Wang [451 146] . as well as Derveaux, Papanicolaou, 
and Tsogka [27]. In [5H [35], an evidence of increasing stability in wave 
imaging when frequency is growing was given. In [15], stability estimates 
for the inverse source problem were established and the conversion of the 
logarithmic type stability to a Lispchitz one first proved. Numerical results 
to illustrate the stability of the source reconstruction problem were pre- 
sented in [16]. In [17], Lipschitz stability estimates for the time-dependent 
wave equation were obtained. In [15], a stability estimate for a linearized 
conductivity problem was derived and its dependence on the depth of the 
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inclusion highlighted. In [36], it is shown that the ill-posedness of the inverse 
acoustic problem decreases when the frequency increases and the stability 
estimate changes from logarithmic type for low frequencies to a Lipschitz 
type for large frequencies. In [27], the enhancement of resolution in the 
near-field was studied and numerically illustrated. Our results in this pa- 
per confirm these important observations in a quite general situation and 
precisely quantify them. 

Our paper is organized as follows. Section Q] is devoted to the stability of 
the reconstruction of the potential from the scattering amplitude (called 
also far-field pattern) in the high frequency regime. Theorem [2] proves that 
the low-frequency component of the potential can be determined in a stable 
way from the scattering amplitude. The threshold frequency determines 
the resolution limit. Section [2] extends the results of Section [1] to the near- 
field measurements. Theorem [3] shows that the same results as those in 
Section [1] hold for reconstructing the potential from measurements of the 
near-field operator. In Section [3] we show that we gain infinite resolution for 
potentials supported near the boundary. If the potential is supported near 
the boundary, then infinite resolution can be achieved in the monotone case. 
Theorem 0] provides a Lipschitz stability result for both the low and high 
frequency components of the potential. Finally, in Appendix El we provide 
useful results on Bessel's functions. 

Finally, we mention that the letter C will be used to denote a universal 
constant which may vary from line to line. We also use A < B to denote an 
estimate of the form A <CB for some constant C. We also use the classical 
notation (x) = y/l + \x\ 2 . 

Acknowledgements. We are very grateful to J. Sjostrand for taking the 
time to discuss some properties of the Laplace transform with us. 

1. Far field pattern 

1.1. Definitions and notations. Let Q be a bounded domain in the Eu- 
clidean space M. d of dimension d > 2, let q G Co°(M rf ) be a real- valued poten- 
tial supported in 0. We use the classical notation D = —id for derivatives 
and consider the Helmholtz equation with potential 

(1.1) D 2 u - \ 2 u + qu = 

at frequency A € Ml := M + \ {0}. Plane waves e 4 ^" 1 ^ propagating along the 
direction oj in S"^ -1 are solutions of the free Helmholtz equation 

(1.2) D 2 u-X 2 u = 0. 

Here denotes the unit sphere in M. d . More generally, plane waves gen- 
erate the set of solutions in the space of tempered distributions, 5'(M d ), of 
the free Helmholtz equation: all solutions of (jl.2p with polynomial growth 
are superpositions of elementary plane waves e l ^ x ' u when uj varies on the 
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sphere S"^ -1 . It will be useful later on to adopt Melrose's notation in [42] to 
designate these solutions: for g G C°°(5' d ~ 1 ), we shall write 

*o(\)g(x) = [ e iX ^g(cj)du. 

Obviously &o(\)g is a solution of (|1.2f) which belongs to iS'(M d ). 

To guarantee the uniqueness of solutions to the Helmholtz equation (jl.ip . 
one can impose conditions on the behavior of solutions at infinity. More 
precisely, we are interested in solutions which can be decomposed 

(1.3) u = + u scat = u m + u scat 

as the sum of an incoming planar wave and a scattered wave satisfying 
Sommerfeld's radiation condition 



(1.4) 



d\x 
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d-1 1 
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->• +CO 




, X 2 / 







uniformly with respect to the direction 9 = A at fixed frequency A £ R?j_. 
The following result from [42 \ Lemma 2.4] holds. 

Proposition 1.1. There exists a unique solution to the Helmholtz equa- 
tion (jl.ip with potential q G C^°(f2) of the form 

(1-5) <p q (x,u,\) = e lX ^ + ^ t (x,io,X), 

where the scattered wave (p q cat satisfies Sommerfeld's radiation condition (jl.4[) 
and which is given by 

(1-6) < at = -R q (\)(e iXx -"q). 

Here R q denotes the meromorphic continuation of the perturbed resolvent. 
Furthermore ip q depends smoothly on (x,uj,\) € R d x S^ 1 x and is 
bounded. 

We choose to denote 



/ ip q (x,uj,\)g(u))duj 
the operator with kernel p q (x,u),\) given by (jl.5p . 

Theorem 1. The scattered wave in the solution (|1.3p to the Helmholtz equa- 
tion (jl.ip given by Proposition 11.11 assumes the form 

e *-Mz| / x \ / i \ 
(1.7) ^( x ) = _^ a J_ W ,A) ^1 as ^+oo, 
s 2 \p| / \|x| 2 / 

where a q is a smooth function on S d ~ l x x and 

^(r'"' A )=-27a(2^) 2 fi<v)<P*<y>"> x )°~**' v *v- 



Proof. We consider the Green function G x corresponding to the free Helmholtz 
equation 

(1.8) (D 2 y -X 2 )G x (x,y) = 8(x-y), 

subject to Sommerfeld's radiation condition (jl.4p . with 5 being the Dirac 
delta function. 

Let R > be large enough so that the ball of radius R contains the support 
of q. By definition of the Green function, for all |x| < R we have 

u{x) = [ (D 2 y -X 2 )G x (x,y)u(y)dy, 
J\y\<R 

and if u is a solution to the Helmholtz equation (11.11) we deduce by Green's 
formula that for any \x\ < R 

u(x) = - G x (x,y)q(y)u(y)dy 

Along the same lines, it is possible to derive a similar identity for the plane 

wave u m (x) = e lXuj ' x 

u in (x) = - J (^(x,y)u^y) - G x (x,y)^^(y)^ da{y), 

taking into account the fact that (D 2 — X 2 )u m = 0. Subtracting the two 
identities gives the following representation formula for the scattered wave 
u scat = u - u in 

u s ™\x) = - G x {x, y) q{y) u{y) dy 

J\y\<R 



\v\=R 



f)G /l 7 ,scat 

-^±(x, y) u scat (y) - G x (x, y)-^ r (y) 



The Green function of the free Helmholtz equation is explicitly given by 



d-2 



(i.9) G ^ = h(-k) 2 fr-vr^HV-Mv-vi), 

where H^j^—i ^ s the Hankel function of first kind and order d/2 — 1 (cf. Ap- 
pendix [X]). The asymptotic behavior of Hankel functions (IA.3j) implies that 
for \y\ large enough 

Since we have 

i i ii x 'y ^n(\ x ^ 2 

\x — y\ = \y\ - -j-r + O' 
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we find for \y\ large enough and fixed A: 



\y\ 

We therefore obtain for fixed x, A and large R = \y\ 

Analogously, we have 

^(x, y) = iXc d (X) lyl-^e^le-^ - (l + O (± 

with c,i(A) = 2^(A/27rz) _ 2~. This leads to 

u scat (x) = - / G x {x, y)q{y)u{y) dy 
J\y\<R 

Sommerfeld's radiation condition implies that the second right-hand side 
term tends to zero when R tends to infinity, so we get 

scat / 



u s ™\ x ) = - J G x (x,y)q(y)u(y)dy. 

Using once again the asymptotic formula (jl.lip together with the fact that 
the Green function is symmetric we get 



To summarize, we have that 

\x\ 2 \\X\ / \\x\ 2 / 



\X\ 2 M-^l / 2 

with 

— ,U),X] = -| I / q(y)u(y)e~ iX ~^' y dy. 

\x\ J 2iX \2iri J J Hyy ' yyj y 

This proves that the scattered part of any solution of Helmholtz' equation 

subject to Sommerfeld's radiation condition takes the form announced in 

Theorem [TJ □ 

Definition 1.2. We define the scattering amplitude associated with the po- 
tential q € C^(R d ) by the smooth function a q : S"^ 1 x S^ 1 xl^C 
given by 

d-l 

(1.12) aq (e,u J ,X) = -^(J^)j 2 jq(y) Vq (y,u J ,X)e^ xe -ydy. 



We denote 

jd-1 



A q (X)g(9) = a q (9,u,X)g{u>)duj, 9eS d ~\ 

Js*- 1 

the corresponding operator with kernel a q (9,u),X). 
It is easy to get an asymptotic expansion of 

%(X)g(x) = I {e iX ^ + ^\x,uj,X))g{uj)duj 

as \x\ — > +00 using the stationary phase and Theorem [TJ 

(1.13) $ q (X)g=(J^j 2 (e-^e^-^gi-e) 



d-l 



with 9 = x/\x\ € S^" 1 . The operator which maps the coefficient of e _ * A ' x ' 
into the coefficient of e lA ' x is given by 



d-l 
2 



(g(0) + (£\ A q {X)g{9)\ . 

This is, after renormalization and composition with the antipodal map, the 
scattering matrix |42j . 



Definition 1.3. The scattering matrix is the operator given by 



d-l 

S q (X) = 1d + (^j 2 A q (X). 



Integration by parts allows to relate the values of the potential inside the 
domain with the scattering matrix S q (X). 

Lemma 1.4. We have the following identities 

J (li ~ <h)u\U2 dx = -2ix(^^j j ^ ^g\~g~2 - S qi (X)g 1 S q2 (X)g 2 dui, 

J (?i - q2)uiu 2 dx = -2iX\^\ J ^ ^ g 2 S qi (X)gi — g\ S q2 (X)g 2 dui, 

for any pair of solution^ 

ui = $ 91 (X)gi , u 2 = $ q2 (X)g 2 , 
to the Helmholtz equations (jl.ip related to the potentials q\,q 2 . 



^We use the notation g((jj) = g(~ijj) = Pg{u) for the antipodal map. 
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Proof. Let R be large enough so that the ball of radius R contains the 
support of both potentials qi,q2- By Green's formula, we have 

/ (<7l — Q2)u\U2 dec = / d r u\ U2 — U\ d r u 2 da(x) 

J\x\<R J\x\=R 

and using the asymptotic formula (|1.13p on the functions u\ = & qi {X)g\ and 
U2 = &q 2 Wd2, we deduce 



d r u\ U2 — ui d r u 2 

d-l 



2iX (jk) (9i{-o)92{-o)-s qi {\) gi {e)s q2 {\)g 2 {e))+o(- 



XRj V yiV "> "WW^wwJJ^^^ndJ' 

which implies 



/ (qi - q 2 )uiu 2 dx 
J\x\<R 



= -2iX ^— 1 J i [g x 52 - S qi9l S q2 g 2 ) d9 + O 

Letting R tend to infinity provides the first identity. 
The proof of the second identity is similar, since 

d r ui u 2 - ui d r u 2 = ~ 2iX {j^) (ai(-o) S q2 g 2 {9) - S qi (\) gi (8) g2 (-e) 

This completes the proof of the lemma. □ 

Choosing q\ = q2 = q in Lemma 11.41 we obtain the following properties of 
the scattering matrix: 

t s q (x) = Pos q (\)op, s q (\ys q (\) = id 

recalling that Pg = g is the antipodal map. Here, t denotes the transpose 
and * the transpose conjugate. Besides, using the first relation together 
with Lemma \1A\ we finally get 

(1.14) 

f 171 — (7oWl ll.n c\t. = OA \ I 

Xi J J S d ~ 



d-l 

(<?i - q2)uiu 2 dx = 2iX( — ) / g 2 (A qi (X) - A q2 (X))giduj, 



for any pair of solutions 

tti = <Z> qi (X)gi, u 2 = $ q2 (X)g 2 , 

to the Helmholtz equations (jl.ip with Sommerfeld's radiation condition (jl.4p 
related to the potentials qx,q 2 . 
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1.2. Stability estimates at high frequencies. The first stability result 
we shall prove in this paper states that low frequencies (i.e. smaller than 2A) 
may be recovered in a stable way from the scattering amplitude. 

Theorem 2. For all e, M, R > and all a > al there exist C £ ,Xo > such 
that the following stability estimate holds true. 

Let q\,q2 £ Cg^R^) ^ e ^ wo potentials supported in the ball centered at and 
of radius R such that ||<72||l°° < M. Then for all A > Ao, 

(i-i5) / (O'lfaT^xe)! 2 ^ 

■'|£|<(2-e)A 

< C £ )?\\a qi — a q2 \\ 2 L 2 + C £ X 2 \\q\ — <?2||i°°- 

The proof of Theorem [2] relies on estimates on solutions to the Helmoltz 
equations stated below, as well as on the following lemma, stating that 
one can relate the Fourier transform of the potential and the scattering 
amplitude. 

Lemma 1.5. For all M,R > there exist constants C, Ao > such that 
for all potentials q G Co°(M d ) supported in the ball B(0,R) and satisfy- 
i n 9 lklU°° — M, we have the following approximation: 



d—B 



< CX^ \\qf L 2 



(1.16) 
for all X > Xq. 

Lemma f 1.5 1 is based on the following classical estimate on the free resolvent 
Rq(X). We refer to [21] for an exposition by Burq of an elementary proof 
due to Zworski. This type of a priori estimates have a long history and go 
back to the work of Agmon [2] on weighted estimates on the resolvent and 
the limiting absorption principle. 

Proposition 1.6. Let x £ Cq°(R ), there exists a constant C > such that 
for all X > 1 we have 

Proof of Lemma \1.5[ We start with the expression defining the scattering 
amplitude (see Definition II. 2h : 

d-1 

a q (0,u,\) = ~(^\ 2 j q(y)<p q (y,uj,X)e- iX9 -ydy 
so by (jl.5p we deduce that 



On 



d-1 

1 / A \~ 



2iX \ 2iri 



q(x)ip s q cat (x,uj,X)e- iXx - e dx. 
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The resolvent identity reads (see for instance Formula (2.3) in 
(1.17) R q (X) = R (X)-R (X)qR q (X) 

which implies that if x € Co°(M d ) is a cutoff function which equals 1 on the 
ball B(0,R), then 

xR g (X)x = xRoWx - {xRo(\)x)q{xR q Wx}- 

If we apply this identity to -qe iXxuj : we get in view of (jl.6p 

X^T = -{xRoWx){qe axuJ ) + {xRoWx) {qX^T) 
and using the estimate of Proposition 11.61 we obtain 

\\X<f%*\\v ^ CX-^qh* +CX~ 1 \\q\\ L ^\\ XV > q cat \\ L 2. 
Taking A > 2CM we deduce 

||x^ cat ||i2<2CA- 1 ||g|| i2 . 
Using Cauchy-Schwarz's inequality, we get 



^sc^-iXx-u dx 



< 2CX- l \\q\\ 2 L 2 



and this completes the proof of the lemma. □ 

Note that the main ingredient in the proof of Lemma [1.51 consists in combin- 
ing the estimate for the non-perturbed resolvent (Proposition II. 6|) with the 
resolvent identity (|1.17|) . along with the fact that q is compactly supported. 

This lemma is the basis of the reconstruction of the potential from the 
scattering amplitude at high frequencies, since one can choose 9, uj € S^" 1 
such that X(6 — u) = £ for any fixed frequency £ with |£| < 2A and let A 
tend to infinity. We use a similar approach to prove a stability estimate at 
high frequencies. 

Proof of Theorem OJ Again, we start with the expression defining the scat- 
tering amplitude and then express the potential in terms of the scattering 
amplitude: for j £ {1,2} we write 

qj(\(0-u)) = -2i{2m)^rX~^a qj (6,uj,X) - j qj(x)(p q x ^ ,t (x, uj, A) e ~ iXx ' e dx. 

Taking the difference of the two expressions yields 

(1.18) \(g^q2){X(e-u))\ <2(2ir)^X-^\(a qi -a q2 )(e,u;,X)\ 

+ 11^. ^,11 II, —SCelt || i 1 1 „ 1 1 ||, _SCclt , _SC3/t 1 1 

lift - <?2||l2 W qi \\L 2 {\x\<R) + llftllL 2 \\<P qi ~ l Pg 2 \\L 2 (\x\<R) 

where we recall that the supports of q\ , q2 are contained in the ball of center 
and radius R > 0. As in the proof of Lemma 11.51 we have 

(1-19) HOIl*(M<*) ^C^Hgilb. 

In light of (|1.6p , the difference of the two scattered waves satisfies 

- <fT = ~(R qi W - R q2 (X))(qie lX ^) - R q2 W((qi ~ ^)e iA ^) 
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and by the resolvent identity (|1.17[) . the first term on the right-hand side 
reads 

fl»(A)(gi " 52)^ 1 (A)(gie iA ^) = -R q2 (X)( gi - q 2 ) ■ 

Let x ^ Co°(]R d ) be a cutoff function which equals 1 on the ball B(0,R), 
then 

x^T - vT) = -(tf?*(A)x) (in - «0<O 

-( x R q2 (\) x )(( qi -q 2 y x ^). 

Therefore along the same lines as in the proof of Lemma 11.5^ the resolvent 
identity (11 . lTf) and Proposition 11.61 give rise to 

ft on\ 1 1,, sca t , .scat II \— 1||„ _ || ||,„scat| 

(1.20) -<p q2 ||i2(|a.|<fl) £ A \\qi - q-2\\L°°\\<P qi \\l 2 {\x\<r) 

+ - Q2\\l 2 - 

Indeed invoking the resolvent identity (|1.17p . we infer that 

xR q2 Wx = xRoWx - {xRo{X)x)q2{xR q2 {X)x) ■ 

Thus applying this identity to (q± — (/2)</?^ at > we obtain 

|| (xR q2 (X)x) (((?i -92)^)11^ < CX- 1 \\q 1 -q 2 \\ Lx y s q ^ t \\ L2(]x] < R) 

+ C\-% 2 \\ L ~\\{ X R q2 (\)x){(qi-q2)vT)\\ L ^ 
which implies for A > 2CM 

\\{ X R q2 (X)x){(qi -92)OHi 2 ^ 2CA- 1 |k 1 - g2 || i ~|Kf ll^(N<«)- 
We have also 

\\( X R q2 (X) X ){(qi-q2)e lXx - UJ )\\ L 2 < A" 1 ^ - q 2 \\ L z, 

which achieves the proof of Estimate (ll.20|) . 

Taking into account those bounds, we obtain in view of f 1 1 . 1 8 1) the estimate 

(1.21) \(q7^2){X(0-uj))\ < \-^\(a qi -a q2 )(9,oj,X)\ 

+ - 92||l°°- 



We denote r(f) = y/1 - |£| 2 when |£| < 1. For f € R d with |f | < 2A, 
choose 7] E with norm \rj\ = r(£/2A). The vectors 

' 2A ' 2A 

have length one, taking the square of (jl.2ip and integrating the resulting 
estimate with respect to 77 yields 

|(gf--?2)(£)| 2 ~ A~ 2 ||gi - g 2 ||ioo 

e e x\ 2 , 

dry. 



+ A-« r «/2A)-« I („„, _„.)(, + ±, ±, A) 
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When |£| < (2 — e)A we have r(£/2A) > y/e(4 — e)/2 therefore multiplying 
the previous inequality by (£) -Q < 1 and integrating with respect to £ we 
get 

/ (CrlteT^XOl 2 ^ < C e A -2 ||qi - gallioc 

■/|f|<(2-e)A 

+ <? £ A 3 y <i j H=r{0 | (o ffl - a, 2 ) (n + £, r? - e, A) | 2 dr, d£. 

We consider the following 2d — 2 dimensional submanifold of S' 2rf_1 

S = {(e,r ? )G5 2 ^ 1 :(£,r / )=0} 
and the following diffeomorphism 

ip : E S^ 1 x 

with inverse 



for which we have 



I F(7 ? + £,7 ? -e)dr ? Ade= / F(0,a;)v9- 1 *(d7 ? AdO 
is Js d - 1 xS d - 1 

= 2~ d [ F{6,uj)duj Ad0. 

Then we finally get 

J|fl<f2-elA 



<C7 £ A 3 // \a gi (e,u,X) -a g2 (e,u,X)\ dOdui 

JJs d - 1 xS d - 1 



+ C e \ ||gi-g2||i«> 
and this completes the proof of our estimate. □ 

In particular, we recover the uniqueness of the potential from the scattering 
amplitude at high frequencies. 

Corollary 1.7. If for all A) belonging to S d ~ l x S d ~ l x we 

have a qi (9,uj, A) = a q2 (6,uj, A) then q\ = q2- 
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2. Near field pattern 



2.1. Definitions and notations. Instead of considering the Helmholtz 
equation on the whole Euclidean space (with Sommerfeld's radiation condi- 
tion) we focus on the Cauchy problem with Robin boundary condition on a 
bounded open set ficl'' with smooth boundary 



This problem has a unique solution u £ H 1 ^) for all / € L 2 (dQ). Writing 
a variational formulation of (|2.ip and using a unique continuation argument 
shows the uniqueness of a solution to (12. 1|) . The existence follows from 
Fredholm's alternative |41j . 

Remark 2.1. Other classical boundary conditions are either Dirichlet or 
Neumann boundary conditions. However, one has to make the additional 
assumption that A 2 is not a Dirichlet (or Neumann) eigenvalue of D 2 + q, 
for the Dirichlet problem corresponding to (12. If) to have a unique solution 
for all / € Hz{d£l). Unfortunately, this does not make sense if one wants to 
take the high frequency limit A —¥ oo. To bypass this difficulty, we study the 
boundary problem (12. ip . Indeed, this condition is natural. It approximates 
Sommerfeld's radiation condition at high frequencies [30|. I38j. 

Remark 2.2. It follows from [H] that the unique solution u £ H l (£l) 
of (12. ip satisfies 



for some constant C independent of A. Therefore, as A — > oo, we have 



The typical inverse problem on the near field pattern is whether it uniquely 
determines the potential q. This was solved (for smooth potentials) by 
Sylvester and Uhlmann |57j in dimension d > 3 in the case of the Dirichlet- 
to-Neumann map. Reconstruction methods were proposed by Nachman 
in |43j and stability issues were studied by Alessandrini [6]. It was shown by 
Mandache in [40] that the logarithmic stability result of Alessandrini in [6] 
is optimal. 



(2.1) 




d u\\l^q) + <MM| L 2(q) < c {\\q u \\l2(q) + ||/||L2(an)) 



DU\\ L 2(Q) + -|M| L 2 (0 .) < C||/|| L 2 (an ). 



Definition 2.3. The near field pattern is the map 



AT q {\) : L 2 (dn) -> L 2 (dn) 

f ^ u\dQ. 



2.2. Stability estimates at high frequencies. The following result is the 
counterpart of Theorem [2] in the near-field context. 
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Theorem 3. For all M, R > and all a > d there exist C, Ao > such that 
the following stability estimate holds true. 

Let qi,q2 € C^°(M d ) be two potentials supported in the ball centered at and 
of radius R such that H^iHl^, H^Hl 00 < M. Then for all A > Ao, 

(2.2) / (e)~ a |gT^2(0| 2 de<C'||A^ 1 -A^ 2 || 2 + C'A- 2 ||( Zl -( Z2 ||| O o. 

J\(,\<2\ 

Proof. We start with two solutions U\, U2 of the equations 

(D 2 - A 2 + qj)uj = 

and computing 

/ (qi ~ ^2)^1^2 dx = / Auiu 2 dx — / uiAu2<ix 
Jn Jn Jn 

{d u — i\)u\U2 da — / u\(d u — i\)u2 da 
an Jan 

yields the formula 

(qi - q 2 )uiU2 dx = [u 2 N qi {X)ui - U\ J\f q2 (X)u 2 ) dcr. 
Jan 

Choosing c/i = c/2 shows that the application J\f q . is symmetric, and using 
this additional information, we get the formula 

(c/i - q 2 )uiu 2 dx = / u 2 (M qi (X) - J\f q2 W)uida. 
Jan 

Let us write q = q\ — q 2 and choose 

ui = <p qi = e iXx ^ + ^ 
and u2 = <p q2 =e iXx - e * + Lp s ™ t , 

with 0i,6»2 € S 11 - 1 . Recall fr om the proof of Lemma 11.51 that the scattered 
waves satisfy the estimate 

1 1 , scat 1 1 ill scat i ^ /~i \ — 1 

W qi \\L 2 {B(0,R)) + Wq 2 \\L 2 {B(Q,R)) S ^A . 

It follows that 

e iXx < e ^q(x)dx= [ (Af qi (\) - M q2 (X)) Ul u 2 da 
n Jan 

„l ™\ f,„scat i ,„scat . scat, _scat\ j_ 
l(X)VP qi + <Pq 2 + l Pq 1 ¥q 2 )dx 

Jn 

and therefore 

C 

\q(\(9i + 2 ))| < \\jVq 1 (X)-jV q2 (X)\\\\u 1 \\ L 2 {dn - ) \\u 2 \\ L 2 {m - ) + - n 

<C\\M qi (X)-M q2 (X)\\ + j\\q\\L^ 
for all (9i,9 2 ) E S^ 1 x S^ 1 . 



V M \\L° 
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We notice that 

S d-i x gd-i _^ 5(0,2) 

(01,02) ^e 1 + e 2 

is a submersion when 9\, 02 are not colinear. This implies that 

C 
A 



(2.3) \q(0\ <C||A^ 1 (A)-A^ 2 (A)|| + -||g||^, £eB(0,2A). 



Multiplying this estimate by (£) a /2, taking the square and integrating on 
the ball -B(0, 2 A) completes the proof of the theorem. □ 

3. The case of a potential located near the boundary 

3.1. Definitions and notations. In this section we show, in the model case 
of the unit disk, that if the potential is supported close to the boundary of 
the disk, then a larger range of frequencies may be recovered by the near 
field J\fq(\) than in the general case treated in the previous two sections. 
More precisely, introducing radial coordinates (xi,X2) = (rcos#,rsin#), 
with (r,6) € M + x [0, 2n[, we consider the following model problem in two 
space dimensions: 



(3.1) 




\ 2 +q\)u n = in B = {x € M. 2 , \x\ < l}, 



We suppose that q\ is a smooth, radial function, with support included 
in D*\ for some fixed constant k > 0, where 

(3.2) Dl = {r £ [0,1], 1 - kA" 1 < r < l}. 

In the following for simplicity we shall drop the index A in the notation of 
the potential. 

The following result shows that in the monotone case, one can improve on 
the frequency band recovered in the general case (see Theorems [2] and [3|) . 
We recall that || JVg(A)|| denotes the operator norm of jV 9 (A) in C(L 2 ). 

Theorem 4. Let q\ and q 2 be two smooth radial potentials supported on D 1 ^ 
as defined in (13. 2ft and such that q\ > q 2 . There are positive constants Xq 
and C such that the following holds. 

Let A i — y be any function such that A < K(\). Then the following 

stability estimate is valid for all A > Aq.' 



i\<K{\) 



| (fi^) (0 1 2 de < CK 2 ( A) (A 4 II M qi ( A) - M q2 (A) |p 



r 2 

, ^qim || i|2 

H ~ " ic 



where C quq2 = max (||<?i||z,°°, ||g 2 |U°° , \\qi \\l<=° \\Q2\\l^ 
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Remark 3.1. One has trivially that 

/ |(gT rr ?2)(0| 2 de<^ 2 (A)||gT^ 2 ||| co 
J\e\<K(\) 

CK 2 (\) U 2 

< ^2 \W - Q^IIl 00 ; 

so the estimate provided in Theorem |4] is of a different nature. 

Remark 3.2. The proof of Theorem U] is presented in Section 13.21 as an 
immediate consequence of Lemma 13.31 proved in Section 13.31 That lemma 
relates the Laplace transform of a function to the near field operator. It holds 
in much more generality than Theorem[4j without the additional assumption 
that q± — q2 > 0. However we are unable to relate the Fourier transform 
of a function to its Laplace transform in general (this fact is well-known to 
be difficult and in general very unstable); for nonnegative functions however 
the relation is very easy and enables us to conclude. 

3.2. Proof of Theorem [4j We start by stating a lemma, proved in Sec- 
tion 13.31 which relates the Laplace transform of a function to the near field 
operator. It is stated in the framework of general, non radial functions. 
We define the operator T : g h-> Tg by 

Tg(r) = g{l - r) 

as well as the Laplace transform £: 

C(g)(t)= [ g(s)e- st ds. 

JR 

For any function 9 i-> f(9) we call Ck(f) its Fourier transform, for k 6 Z: 

c k (f)(r)= f(r,9)e-* kd d6. 
Jo 

Lemma 3.3. Let q\ and q<i be two functions supported on D 1 ^. There are 
three positive constants \q-,K and C such that ift> 2KX > 2KXq, then 



£{c k {T(r( qi -q 2 )))(t) 

< c (\ 2 \\M qi (X) -MMW + ^Iki - <? 2 ||l° 



A 2 

where C qim = max (||gi||z,°°, IMIl^, ||gi||L°°|M|i,°°)- 

Proof of Theorem^ Let q = q\ — q 2 - Since q is radial, we have q(x) = Q{\x\) 
and 

m) = Q(iei) 

with for all p > 0, 

r 2ir roc 



rzn /'OO 

Q(p)= / e- irpcosd Q(r)rdrd9. 
Jo Jo 



In particular recalling that Q is nonnegative, we find that for all £ € 

/>oo 

< 2tt / Q(r)rdr 
J o 

/>oo 

<2tt / e ?o(1 - r) e- Co{1 - r) Q(r)rdr 
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p2 



< Ce 3KK ( A 2 \\N qi {\) -A^ 2 (A)|| + -^||gi - <7 2 ||l~ 



for any (o € R. Then we can apply Lemma 13.31 to the particular case of a 
radial function, so choosing k = and Co = 3-ftTAo we infer that 

(3.3)|g(OI < Ce^ K / A £(T(rQ))(Co) 

A 2 

The end of the proof of Theorem [J] follows easily by taking the L 2 norm 
in £. □ 

3.3. Proof of Lemma 13.31 The method of proof follows the ideas devel- 
oped in the proofs of Theorems [2] and [SJ adapting the estimates to our special 
situation where the potentials are located near the boundary. The heart of 
the matter consists to approximate the solutions to Helmholtz equation ([3.1 1) 
by separable solutions in radial coordinates involving Bessel functions, which 
allows in light of Debye's formula to relate the Laplace transform of the po- 
tential to its near field operator. 



More precisely, we shall look for solutions to ([3JJ) under the following 
form, for n € Z: 

/ a\ J |n|(Ar) . Q 

X ( J \n\( X ) - tJ \n\W) 

= z n (r, A) e ind + v n , 
where J n is a Bessel function of the first kind (see Appendix [A}, solution to 

Suppose that for I £ {1, 2}, u^ n solves f)3. 1 j) with potential qi. Then writing 

u e n (r,e)=z n (r,X)e me + v e n (r,d), 

we find that 

((D 2 -X 2 + q e )vi = -q e z n e ine in B, 

Due to Property (|A.6|) we have if 1 — A" 1 < r < 1 and \n\ > KX for K large, 
that \Xz n (r, A)| is bounded, for A > Ao, by a constant depending only on K 
and Ao- 

For now on we shall denote by C(Ao, K) such a constant, which may change 
from line to line. 
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Therefore, arguing as in the proof of Lemma 11.51 and taking advantage of 
the fact that on the support of qt, r varies in an interval of size A -1 , we find 
as soon as A is large enough compared to ||<#||l°° 



(3.4) 



\v n \\ L 2 < \\q e {r,e)z n {r,X)\\ L 2 < 



A 



A§ 



Now going back to the computations of Section [2] we consider two positive 
integers n and m such that n, m > KX, and we write 



( qi - q 2 )Kut m dx = / (M qi (A) - M q2 (A)e my e" imy d9. 

B JdB 



Therefore decomposing 



uKr, 9) = z n (r, A) e md + t£(r, 9) and 
u 2 _ m (r, 9) = z m (r, A) e~ imd + v 2 _ m (r, 9), 



we get 



2tt 



>, 1 (A)-AT (?2 (A)) 



(qi ~ g 2 )e i(n - m)9 z„(r, A) z m (r, A) rdr d9 



(3.5) 



B 



+ / (qi-q2)e in8 z n (r,X)v 2 _ m rdrd9 
Jb 

+ [ (qi-q2)e- ime z m (r,X)v 1 n rdrd9 
Jb 

+ / (<?i - Q2)vnV 2 . m rdrd9. 

JB 



As Az n (r, A) | is bounded by a constant C(Aq, K), we have (recalling that q\ 
and q2 are compactly supported in an interval in r of size A -1 ) 



(qi-q2)e iny z n (r,X)vi m rdrd9 
and similarly 

(qi - q2)e- ime z m (r, X)v 1 n rdrd9 



IB 
Finally 



C(X ,K) 2 

X2 



< — i — hi - q2\\L°°\\v n \\ L 2. 

A2 



(qi - q 2 )v n v_ m rdrd9 



< hi - 92||l°°II^||l2||^ 



Ii2. 
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Thus by MM . we get from (pO)|h for n,m> KX> KX , 



(3.6) 



(ffi - (? 2 )e i(ri - m)e z n (r, A) z m (r, A) rdr dfl 



< 



2tt 



> 9l (A)-^ g2 (A)) 



e ine e -ime de 



+ 



G 



A 4 



q\\L° 



where q = q\ — qi and 

£91,92 = C(X ,K)max (||gi||i<x», ||g 2 |U°° , ||?i||z,°° IM|l°°)- 
This gives rise to 



(3.7) 



I (?l ~ q2y (n - m)e z n (r, X) z m (r, A) rdr d0 
Jb 



< WfV qi (\)-fV q2 (\)\\ + 



a 



A 4 



q\\l°°, 



which leads in light of (|A.6j) to 

(qi ~ q 2 )e i{n - m)9 e- in+m){1 - r) rdrdO 



B 



< C(A ,K)(a 2 ||AUA) - AUA)|| + 



91.92 



A 2 



\<1\\l° 



for all n, m > K X. In conclusion we have for any k 6 Z, any £ € N and 
any j > 2KA > 2KA , 



(3.8) 



/ (91 -g 2 )e- ifc V j(1 - r) rdrdfl 
JB 



91.92 



<C(A ,K)(a 2 ||A^ 1 (A)-AA (?2 (A)|| + 
The conclusion follows from (|3.3p . Lemma 13.31 is proved. 



□ 



4. Concluding remarks 



In this paper we have shown that the low-frequency component of the 
potential can be determined in a stable way from the scattering measure- 
ments and justified the resolution limit. We have also proved that in the 
near-field we have in the monotone case infinite resolution in reconstructing 
the potential near the boundary. We think that the result holds in the gen- 
eral case. However, its proof seems to be out of reach. In fact, even though 
a sampling (or interpolation) formula for the Laplace transform does ex- 
ist |20l [54] , making norm-estimates similar to those in Theorem |4] is very 
challenging. Our results can be extended in many directions. It would be 
very interesting to study the limited- view case and show, as in [IT], that 
we recover infinite resolution from near-field measurements on the overlap 
of the source and receiver apertures. Another challenging problem is to un- 
derstand how probe interaction can improve local resolution by converting 
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evanescent modes of the potential to propagating ones [23j . These problems 
will be the subject of forthcoming works. 

Appendix A. Bessel functions 

Bessel's equation arises when finding separable solutions to the Helmholtz 
equation in spherical coordinates, and writes as follows: 

(A.l) (d 2 z + -d z + (l-\))u = 0, neZ. 



z 



It is well known (see for instance [391 09] and the references therein) that one 
of the solutions of Bessel's equation is the entire function J n (z) known as the 
Bessel function of the first kind of order n, and defined for arbitrary z £ C 
by the convergent series 

k=oo (_-l\k 

k=0 y ' 
in the case where n € N and by J- n (z) = (— l) n J n (z). 

To find a general solution of Bessel's equation (jA.ip . we need a second 
solution of (jA.ip which is linearly independent of J n (z). For such a solution, 
we usually choose Y n (z) the Bessel function of the second kind which is entire 
in the complex plane cut along the segment ] — oo, 0] and defined for n € N 

by 

7T 2 7T k\ 2 

k=0 



k=oo 



7r ^ k\(n + k)\ 

where ib(m + 1) = —7 + 1 H h ... H , 7 being the Euler constant. We 

2 m 
also define Y_ n (z) = {-l) n Y n (z). 

Since J n and Y n are linearly independent, the general expression for solutions 
of (|A.1|) is a linear combination of Bessel functions of the first and second 
kinds, i.e, 

u{z) = A J n (z) + B Y n (z), 
where A and B are constants. 

Another basis of solutions to the differential equation (jA.lh is given by the 

Bessel functions of the third kind or Hankel functions, denoted by H^p 

(2) 

and Hn . These functions are defined by the formulas 

(A.2) H^(z) = J n (z)+iY n (z) and H%\z) = J n (z) - iY n {z), 
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where z is any point of the complex plane cut along the segment ] — oo, 0]. 
The motivation for introducing the Hankel functions is that the linear com- 
bination of J n (z) and Y n (z) have very simple asymptotic expansions for 
large \z\: it is thus well-known that 

(A.3) HiP(z) = (J-)V<*-¥-f> (l + O^)) and 
(A.4) . (A)V.(^-f)( 1 + (^)) 

as \z\ — > 00. 

Furthermore, we have the following Debye formulas whose proof can be 
found for instance in |49t Chapter 9.4] and [TJ Chapter 9]: 

g— n(a— tanh a) / ^ g 

J n (n secha) ~ - I 1 + (— cotha (cotha) 3 ) — h 

(27m tanh a) 2 V 8 24 n 

and particulary 

g— n(a— tanh a) / , j 

J n (?i secha) = r 1 + 0[- 

(27rntanha)2 V ^ n 

(A 5) 

, , (tanha)^e-"( a - tanha ) / 
4(n secha) = ^ - ; 1 + O 1 

(47rn)2 V 



as n — > 00, where sechz denotes the hyperbolic secant of z defined by 

sechz 



cosh z 

with coshz the hyperbolic cosine. 

Debye's formula gives rise to the following asymptotic behavior for the func- 
tion introduced in Section 

*(A,r)- J "" (Ar) 



MJ' M (\) - u lnl (Xj) 

defined for n £ Z. Let us prove that for |n| S> A and r ~ 1, 



(A.6) 



and Im z n (A, r) - ^ e - n(1 " r) (l + ofj^j. 



Without loss of generality we may assume that n £ N, then defining 

71 Tl 

coshai = — and cosh = — , 
Ar A 
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it is easy to see that under the above assumptions (A <C n and r ~ 1), 
we have necessarily cosher 3> 1 for i £ {1,2}, hence «j > 1. This im- 
plies cosh a i ~ e°% sinhaj ~ e Qi , and tanhaj ~ 1, which gives rise to 

J n (n sechai )J' n {n secha^) 



Rez n (X,r) 
lmz n (\,r) 



A(J,j 2 (n secha^) + Jn( n secha2)) 

and 

J n (n sechai )J n (n secha^) 



X(Jn(n sechai) + Jn{n secha^)) 
Finally, taking advantage of (|A.5P we get 



Re z n (X, r) 



e -n(a 1 -a 2 ) g-nlog(l) 



A A 

The computation is identical for Imz n (A,r), so using the fact that r is near 
to 1, we obtain the desired conclusion. 
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